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Magneto-solid-mechanical model of two-component, core-crust, paramagnetic neutron star re- 
sponding to quake-induced perturbation by differentially rotational, torsional, oscillations of crustal 
electron-nuclear solid-state plasma about axis of magnetic field frozen in the immobile paramag- 
netic core is developed. Particular attention is given to the node-free torsional crust-against-core 
vibrations under combined action of Lorentz magnetic and Hooke's elastic forces; the damping is 
attributed to Newtonian force of shear viscose stresses in crustal solid-state plasma. The spectral for- 
mulae for the frequency and lifetime of this toroidal mode are derived in analytic form and discussed 
in the context of quasi-periodic oscillations of the X-ray outburst flux from quaking magnetars. The 
application of obtained theoretical spectra to modal analysis of available data on frequencies of os- 
cillating outburst emission suggests that detected variability is the manifestation of crustal Alfven's 
seismic vibrations restored by Lorentz force of magnetic field stresses. 
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I. INTRODUCTION 

The investigations of neutron star seismic vibrations 
offer unique opportunity of studying their internal struc- 
ture, solid-mechanical and electrodynamical properties 
of superdense degenerate matter. The most conspicu- 
ous feature of these non-convective solid stars is the ca- 
pability of accommodating magnetic fields of extremely 
high intensity [l| that serve as a chief promoter of their 
observable electromagnetic activity. The absence of nu- 
clear energy sources in these final stage (FS) stars sug- 
gests that their magnetic fields are definitely not gener- 
ated by persistent current-carrying flows in self-exciting 
dynamo processes, as is the case of liquid main-sequence 
(MS) stars. It seems quite likely, therefore, that stability 
to spontaneous decayQ of fossil magnetic fields of iso- 
lated neutron stars [3| is maintained by permanent mag- 
netization of neutron-dominated (poorly conducting) de- 
generate Fermi-matter. Such an understanding has been 
laid at the base of paramagnetic neutron star model 0- 
7]. In this model the degenerate Fermi- matter of non- 
relativistic neutrons (whose degeneracy pressure with- 
stands the pressure of self-gravity) is regarded as being in 



the permanently magnetized state of field-induced Pauli's 
paramagnetic saturation which is characterized by align- 
ment of spin magnetic moments of neutrons along the 
axis of frozen-in magnetic field. The most striking dy- 
namical manifestation of spin paramagnetic polarization 
of non-conducing neutron matter is that such a mat- 
ter can transmit perturbations by transverse magneto- 
mechanical waves; in such a wave the vector-fields of 
magnetization and material displacements undergo cou- 
pled differentially rotational vibrations traveling along 
the axis of magnetic field. In a spherical mass of param- 
agnetic neutron star this unique feature of field-induced 
spin magnetic polarization of neutron matter is mani- 
fested in that such a star can undergo solely torsional 
vibrations about axis of its dipole magnetic moment. 
Based on this finding, it was argued in above works 
that the model of paramagnetic neutron star execut- 
ing torsional axisymmetric vibrations, weakly damped 
by nuclear matter viscosity, is able to explain long pe- 
riodic ([5 < P < 12] s - non- typical to young neutron 
stars) pulsed character of magnetar radiation (both, Soft 
Gamma Repeaters and Anomalous X-ray Pulsars) in seis- 
mically quiescent regime of their emission, as being pro- 
duced by torsional vibrations, rather than rotation as is 
the case of radio pulsars. 
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Recent years have seen a resurgence of interest 



in torsional vibrations of magnetars, prompted by 
observations [§-[To| of quasi-periodic oscillations (QPOs) 
during the outburst flare from SGR 0525-66, SGR 1806- 
20, and SGR 1900+14. The statistics of X-ray burst of 
SGRs exhibits typical for earthquakes features 11]. It is 
believed, therefore, that detected QPOs are of seismic 
origin. Particular attention in this development of mag- 
netar asteroseismology has been paid to the following set 
of data on QPO frequencies fl2| 

SGR 1806 - 20 : 18, 26, 29, 92, 150, 625, 1840; (1) 
SGR 1900+ 14 : 28,54,84,155 [Hz]. (2) 

The corresponding periods are substantially shorter than 
the above mentioned periods of seismically quiescent 
pulsed emission. In works [l3l - [T7| motivated by this dis- 
covery, several models of post-quake vibrational relax- 
ation of above magnetars have been investigated. Partic- 
ular attention has been given to the regime of node-free 
or nodeless shear axisymmetric vibrations. This regime 
is interesting in its own right because such vibrations 
have been and still are poorly investigated in theoretical 
asteroseismology of both solid FS stars and such solid ce- 
lestial ob jects as Earth-like planets (l^.[l9j. In particular, 
in works [13l - [l5l |. a case of the elastic- force-driven node- 
less shear oscillations, both torsional - otg and spheroidal 
- qS£ , entrapped in the crust of finite depth AR has been 
studied in some details with remarkable inference that 
dipole overtones of spheroidal and torsion vibrations of 
crust against immobile core exhibit features generic to 
Goldstone soft modes. On the other hand, one can prob- 
ably cast doubt on arguments of the model presuming the 
dominant role of solid-mechanical Hooke's force of elastic 
stresses because such interpretation rests on poorly jus- 
tifiable assumption about dynamically passive role of an 
ultra strong magnetic field, that is, that the field frozen 
in the star remains unaltered in the process of vibrations. 
Bearing this in mind and assuming that the presence 
of charged particles in neutron-dominated stellar mat- 
ter imparts to it the properties of electric conductor, in 
works the post-quake relaxation of above magne- 

tars has been studied in the model of perfectly conducting 
solid star executing global torsional vibrations restored 
by joint action of Lorentz force of magnetic field stresses 
and Hooke's force of solid-mechanical elastic stresses. 
It was found that such a model provides fairly reason- 
able account of general trends in QPO frequencies for 
all data from SGR 1900+14 and for SGR 1900+14 from 
the range 30 < v < 200 Hz, but faces serious difficul- 
ties in interpreting low-frequency vibrations with v = 18 
and v = 26 Hz in data from SGR 1806-20. Also, the 
model of global torsional vibrations leaves some uncer- 
tainties regarding the nature of vibrations with v = 625 
and v — 1840 Hz. This last issue has been scrutinized 
in recent work pj] from the standpoint of a solid star 
model with non-homogeneous poloidal magnetic field of 
well-known Ferraro's form. And it was found that these 
high-frequency QPOs can be properly explained as being 
produced by very high overtones of node-free torsional 



Alfven oscillations. 

PARAMAGNETIC NEUTRON STAR 




FIG. 1. (Color online) The internal constitution of two- 
component, core-crust, model of paramagnetic neutron star. 
The massive core is considered as a poorly conducting per- 
manent magnet composed of degenerate Fermi-gas of non- 
relativistic neutrons in the state of Pauli's paramagnetic satu- 
ration caused by field-induced alignment of spin magnetic mo- 
ment of neutrons along the axis of uniform internal and dipo- 
lar external magnetic field frozen in the star on the stage of 
gravitational collapse of its MS progenitor. A highly conduct- 
ing metal-like material of the neutron star crust, composed 
of nuclei embedded in the super dense degenerate Fermi-gas 
of relativistic electrons, is regarded as electron-nuclear solid- 
state magneto-active plasma capable of sustaining Alfven os- 
cillations. 



As a logical extension of above line of investigation, in 
this paper we consider in some details a case of node- 
free torsional vibrations locked in the crust with focus 
on toroidal Alfven mode. In so doing we work from the 
two-component model of paramagnetic neutron star, pic- 
tured in Fig.l, whose crust and core materials are re- 
garded as endowed with substantially different electro- 
dynamic properties. The immobile massive core, pri- 
marily consisting of degenerate neutron matter in the 
above described permanently magnetized state of Pauli's 
paramagnetic saturation, is regarded as a main source 
of magnetic field of the star crust. This implies that 
the core material is just incapable of sustaining Alfven 
vibrations which owe their existence to extremely large 
(effectively infinite) electrical conductivity of matter [20l— 
|22j. The micro-composition of crust, which is dominated 
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by nuclei embedded in degenerate Fermi- gas of rela- 
tivistic electrons, suggests that its metal-like material 
possesses properties of perfectly conducting solid-state 
plasma. Such a view suggests that seismic stability of 
the star to quake-induced tectonic displacements of crust 
against core is primarily determined by well-known effect 
of magnetic (magnet-metal) cohesion mediated by mag- 
netic field lines which operate as a super-hard piles en- 
dowing the core-crust construction of neutron star with 
supplementary (to gravity forces) stiffness of magnetic 
nature. The intermediate layer between core and crust 
(the inner crust whose density several times less than the 
core density) is most likely composed of quasi-boson mat- 
ter of paired neutrons. But it is highly unlikely that such 
quasi- boson matter is capable of undergoing phase tran- 
sition to the Bose-Einstein Condensation (BEC) which 
is characterized by vanishingly small pressure. This sug- 
gests that BEC state of paired neutrons, if exist, can 
only insufficiently contribute to the total mass budget of 
neutron star - a compact object in which pressure of self- 
gravity is brought to equilibrium by degeneracy Fermi- 
pressure of non-relativistic neutrons in the core and rela- 
tivistic electrons in the crust. In seismo-dynamics of the 
paramagnetic neutron star under consideration the inner 
crust is thought of as operating like a lubricant facilitat- 
ing differentially rotational shear displacements of crust 
relative to much denser matter of massive core. From 
the view point of this core-crust model, the star-quake is 
though of as impulsive release of energy of magnetic core- 
crust cohesion (by means of disruption of magnetic field 
lines on the core-crust interface) resulting in the crust 
fracturing by revealed magnetic stresses. In this paper 
we focus, however, not on dynamics of quake, but on 
the post-quake vibrational relaxation of the star, namely, 
on node-free torsional oscillations of crustal solid-state 
plasma about axis of magnetic field frozen in an immo- 
bile paramagnetic core. In section 2, a brief outline is 
given of theory of solid-magnetics appropriate for the 
perfectly conducting viscoelastic continuous medium per- 
vaded by a magnetic field. In section 3, the spectral 
formulas for the frequency and lifetime of differentially 
rotational, torsional, nodeless vibrations of the crust re- 
stored by combined action of magnetic Lorentz and elas- 
tic Hooke's forces are obtained. In section 4, the com- 
puted frequency spectra are used for the forward astero- 
seismic analysis of the fast oscillations of X-ray outburst 
from above mentioned magnetars. The obtained results 
are briefly summarized in section 5. 



deal with the shear differentially rotational fluctuations 
of viscoelastic crustal matter of density p which are de- 
scribed by quake-induced material displacements it; (ba- 
sic variable of solid-mechanics). The non-compressional 
character of vibrations under consideration implies that 
Sp = — pVfeUfc = 0. With this in mind, the govern- 
ing equations of solid-magnetics (solid-mechanical coun- 
terpart of equations of magneto-fluid-mechanics) can be 
written in the form 

pui = VkT ik + VfcO^fe + V k ^ik, V fc u fe = (3) 

presuming that Hooke's elastic stresses and Newton's 
viscous stresses 7r^ are described by linear constitutive 
equations 



Oik — 2p,Uik, u ik 
TT ik = 2rju lk 



1 



[VjUfc + VfcitJ, 



u ik = -jpJiUk + Vfciti 



(4) 
(5) 



where p stands for the shear modulus, rj for shear viscos- 
ity and Uik is the tensor of shear strains or deformations. 
The central to our further discussion is the tensor of fluc- 
tuating magnetic field stresses 

m = ^- [BiSB k + B k 5Bi - BjSBjSik] , (6) 

47T 



SBi = V fe [ui B k - u k Bi] 



(7) 



As in our previous works[16ll29j. we consider model with 
homogeneous internal magnetic field whose components 
in spherical polar coordinates read 



B r = B cos ( 



B e = -B sin f 



B^ = (8) 



and external dipolar magnetic field is described by B = 
V x A, where A = [0,0, A$ — m s /r 2 ] is the vector po- 
tential with the standard parametrization of the dipole 
magnetic moment m s = (1/2)BB? of star of radius R 
and by B is understood the magnetic field intensity at 
its magnetic poles, B = B p . 



A. The energy method 

This method of computing frequency of shear vibra- 
tions rests on the equation of energy balance 



at 



pu 



dV 



Jik + O lk + TT ik } u ik dV (9) 



II. GOVERNING EQUATIONS OF 
SOLID-MAGNETICS 



It is generally realized today that seismic vibrations 
of superdense matter of non-convective FS-stars (white 
dwarfs, pulsars and quark stars) can be properly de- 
scribed by equations of solid-mechanical theory of vis- 
coelastic continuous media |23h28I. In what follows we 



which is obtained by scalar multiplication of equation 
of magneto-solid-mechanics ([3]) with Ui and integration 
over the volume of seismogenic layer. From the technical 
point of view, the shear character of material distortions 
brought about by forces under consideration owes its ori- 
gin to the symmetric form of stress-tensors in terms of 
which these forces are expressed. It is this last feature 
of solid-mechanical elastic stresses and magnetic field 
stresses that endows the solid-state plasma pervaded by 
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homogeneous magnetic field with the capability of re- 
sponding to non-compression perturbation by reversal 
shear vibrations (which are not accompanied by fluctua- 
tions in density). The physical significance and practical 
usefulness of the energy method under consideration is 
that it can be efficiently utilized not only in the study of 
non-radial seismic vibrations of neutron stars, but also 
can be applied to the study of more wide class of solid 
degenerate stars like white dwarfs stars [13, HH and ul- 
tra dense quark-matter stars 32J whose material is most 
likely in the solid aggregate state [H, HIl- At this point it 
seems appropriate to mention here theoretical investiga- 
tions of vibration properties of atomic nuclei (thought of 
as ultra fine pieces of continuous nuclear matter) in which 
it has been found that nuclear giant-resonant excitations 
(fundamental vibration modes generic to all nuclei of pe- 
riodic chart) are properly described in terms of spheroidal 
and torsional elastic vibrations of a solid sphere 35] . This 
suggests that degenerate nucleon Fermi-matter, regarded 
as continuous medium, can be thought of as a strained 
Fermi-solid, rather than flowing Fermi-liquid. 

The key idea of this method consists in using of the 
following separable form of material displacements 



Ui(r,t) = ai{r)a(t) 



(10) 



where Oj(r) is the time-independent solenoidal field and 
amplitude a(t) carries information about temporal evo- 
lution of fluctuations. Thanks to this form of Uj, all the 
above tensors of fluctuating stresses and strains take sim- 
ilar separable form 



T ik (r,t) = [fi fe (r) - -fjj(r)Sik]a(t), 



(11) 



f tt (r) = —[Bi(T)b k (T) + B k (T)bi(T)], (12) 

bi(r) = V fc [oi(r) S fe (r) - a fe (r) B t (r)}, (13) 
cr ife (r) = 2fia ik (r)a(t), ir ik (r) = 2i] a ik (r)a(t), (14) 



«ifc( r ) = 2 [ V ^( r ) + VfcOi(r)]. 



(15) 



On inserting (|TU| - p3)) in the integral equation of energy 
balance ([5]) we arrive at equation for a(t) having the well- 
familiar form 



d£ „^ „ Ma 2 JCa 2 ^ Va 2 

^ = - 2T > £ = — + — T =— (16) 



Ma + Va + Ka = 0, 

a(t) = cto cxp(— t/r) cos(fit), 



(17) 
(18) 



K 



2M 



n 2 =lo 2 [1- [urY 2 ] ," 2 = jj, t=—. (19) 



V 



of damped oscillator are given by 

M = J p(r)ai(r) <Xj(r) dV, K = K e + K, m , (20) 

K e = 2 J fi(r)a ik (r)a ik (r)dV, (21) 

Km= I T lk {r)a lk (r)dV, (22) 



V 



rj(r) a ik (r) a ik (r) dV. 



(23) 



where the inertia Ai, viscous friction T> and stiffness K, 



All the above equations valid for arbitrary volume occu- 
pied by magneto-active solid-state plasma whose density, 
shear modulus and shear viscosity are arbitrary functions 
of position. As in our previous works, here we confine 
our computations to the case of uniform profile of these 
later parameters. And because the purpose of our present 
study is the frequency spectrum of node-free torsion vi- 
brations about magnetic axis of the star, several com- 
ments should be made regarding the axisymmetric field 
of material displacements a which is taken in one and 
the same shape in computing parameter of inertia M 
and spring constants of both solid-mechanical JC e and 
magneto-mechanical IC m stiffness. 

From the physical point of view, the main argument 
justifying the use of one and the same field of mate- 
rial displacements in computing frequency spectra of tor- 
sional vibrations driven by forces of elastic and magnetic 
field stresses rests on general statement of continuum- 
mechanical theories of magneto-active perfectly con- 
ducting continuous media - magnetic field pervading 
(both liquid-state and solid-state) plasmas imparts to 
such a medium a supplementary portion of elasticity 
which is manifested in its capability of transmitting non- 
compressional mechanical perturbations by transverse 
Alfven waves. Such a view is substantiated by the com- 
monly known fact that transverse wave in incompressible 
continuous medium is the feature of material oscillatory 
behavior which is generic to elastic solid, not an incom- 
pressible flowing liquid. The transverse hydromagnetic 
wave propagating, along the lines of constant magnetic 
field B frozen-in the perfectly conducting medium, with 
Alfven speed v A = [2P B /p] 1/2 (where P B = B 2 /8tt is the 
magnetic field pressure), is characterized by dispersion 
equation u> — va k which is similar to that for transverse 
wave of shear mechanical displacements, oj = Cj k, travel- 
ing in an elastic solid with the speed Ct = [n/p] 1 ^ 2 where 
p is the shear modulus (which has physical dimension of 
pressure). The profound discussion of analogy between 
oscillatory behavior of incompressible perfectly conduct- 
ing plasmas and elastic solid (similarity of transverse hy- 
dromagnetic wave in incompressible perfectly conducting 
media and transverse wave of shear mechanical displace- 
ments in an elastic solid) can be found in monographs of 
Chandrasekharf^ (section Alfven waves), and more ex- 
tensively this issue is discussed in monograph of Alfven 
and Faltammar[21]. Regarding the difference between 
node-free oscillatory behavior of solid sphere and spher- 
ical mass of an incompressible liquid it is appropriate 
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to note that the liquid sphere is able to sustain solely 
spheroidal node-free vibrations of fluid velocity. The 
canonical example is the Kelvin fundamental mode of 
oscillating fluid velocity in a heavy spherical mass of in- 
compressible homogeneous liquid restored by forces rep- 
resented as gradient of pressure and gradient of potential 
of self-gravity [3^| . In the meantime, the node- free vibra- 
tions of solid sphere restored by elastic force (represented 
as divergence of shear mechanical stresses) are charac- 
terized by two eigenmodes. Namely, the even-parity 
spheroidal mode of non-rotational vibrations of material 
displacements and the odd-parity torsional mode of dif- 
ferentially rotational vibrations, the problem in which the 
very notion of the torsion vibration mode has come into 
existence fl3l j . In our studies focus is laid on poorly inves- 
tigated regime of node-free vibrations in which solenoidal 
field of material displacements, V • a = 0, obeys the vec- 
tor Laplace equation, V 2 a — 0. In this regime the in- 
stantaneous material displacements are described by the 
toroidal field of the form a = AgV x [rr Pi (cos 6*)]. Sub- 
stituting this field in the above given integrals for inertia 
Ai and solid-mechanical stiffness /C e and integrating over 
the entire volume of oscillating star we have found in our 
previous studies that elastic-force-driven node-free global 
torsional vibrations, M.a+K. e a = 0, are characterized by 
the frequency spectrum ui e (£) — [K, e / M] 1 / 2 of the form 

u 2 (t)=u 2 l(2e + 3)(l-l)], (24) 



UJ,, 



R' 



Ct 



(25) 



where uj e is the natural unit of frequency of shear elastic 
vibrations and c t is the speed of transverse wave in the 
elastic solid characterized by average shear modulus /2 
and average density p. 

The above line of argument about similarity between 
oscillatory behavior of elastic solid and magneto-active 
plasma suggests that perfectly conducting matter of neu- 
tron star pervaded by homogeneous magnetic field should 
be able to sustain the Lorentz-force-driven differentially 
rotational seismic vibrations (triggered by quake) about 
magnetic axis in which oscillating field of material dis- 
placements has one and the same form as in the above 
outlined Hooke's-force-driven torsional vibrations. Ad- 
hering to this assumption and making use of the above 
node-free toroidal field, a, as a trial function for comput- 
ing A4 and K, m we found in [HI [ItJ that torsional vibra- 
tions restored by magnetic Lorentz force, Aia+JC m a = 0, 
are characterized by the frequency spectrum 
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VA 

R' 



B 



VA 



(26) 
(27) 



where ua is the natural unit of frequency of Alfven vi- 
brations and va is the speed of Alfven wave. The prac- 
tical usefulness of outlined computations with one and 
the same trial toroidal field of displacements is that al- 
lows us to assess the relative role of restoring Hooke's and 



Lorentz forces in torsional seismic vibrations of neutrons 
stars which are responsible, as is believed, for the fast os- 
cillations of X-ray flares from quaking magnetars. With 
all above in mind, one of the main purposes of present 
paper is to make such an assessment in a mathematically 
consistent fashion for the torsional node-free vibrations 
entrapped in the crust, the problem which is considered, 
to the best of our knowledge, for the first time. 



B. Material displacements in torsional mode of 
crust-against-core nodeless vibrations 

The above equations of the energy method show that 
the main trial function of the frequency spectrum com- 
putation is the toroidal field of instantaneous displace- 
ments. Because one of the main our purposes here is to 
assess the relative role of elastic and magnetic forces in 
quake-induced torsion nodeless vibrations of magnetars, 
we again adopt all the above arguments regarding the 
choice of this field in the form of the general solution to 
vector Laplace equation. 

It is convenient to start with the rate of material dis- 
placements which is described by general formula of ro- 
tational motions 

Hr,i)=u(r,t) = [n( r) t)xr], (28) 
f2(r,t) = [V x 5v(r,t)} = [V x u(r,t)] = *(r,*)(29) 

However, unlike a case of rigid-body rotation, in which 
the angular velocity is a constant vector, in a solid mass 
undergoing axisymmetric differentially rotational vibra- 
tions the angular velocity ft(r,t) is the vector-function 
of position which can be represented as 

n(r, t) = *(r, t) = 4>{r) a(t), 0(r) = [V x a(r)] (30) 

In the regime of node-free vibrations in question, a(r) 
is described by the divergence-free odd-parity, axial, 
toroidal field which is one of two harmonic solenoidal 
fields of fundamental basis [2^ obeying the vector Laplace 
equation V 2 a = 0. This field can be expressed in terms of 
general solution of the scalar Laplace equation as follows 

a(r) = a t (r) = V x [r X (r)] = [V X (r) X r] (31) 
V 2 x(r) = 0, (32) 
X (r) = [A t r l + B t r- l - 1 ]Pt{Q, ( = cos0 (33) 

where Pt(C) is the Legendre polynomial of multipole de- 
gree £. It follows that the angular field 0(r) is the 
poloidal vector field 

</>(r) = [V x a t (r)] = V x V x [rx(r)] (34) 
= V[A e r e + Btr- e - 1 ]P t (Q, (35) 
At=Ai{l + l), B f = B e £. (36) 

And this field is irrotational, V x <p(r) = 0. 

As was stated, we study a model of differentially rota- 
tional vibrations of peripheral finite-depth crust against 
immobile core. In this case, the arbitrary constants Ae 
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FIG. 2. (Color online) Material displacements of crustal mat- 
ter about the dipole magnetic moment axis of paramagnetic 
neutron star undergoing nodeless differentially rotational, tor- 
sional, vibrations in quadrupole and octupole overtones. 



and Be can be uniquely eliminated from two boundary 
conditions: (i) on the core-crust interface r = R c 



U<p\ r =R c = 

and (ii) on the star surface r = R 

U<j,\ r =R = [* X R]0| r =R, 



where 



1 d 



oO sin ( 



(37) 



(38) 
(39) 



(40) 



The no-slip condition on the core-crust interface, r = R c , 
reflects the fact that the amplitude of differentially rota- 
tional oscillations gradually decreases down to the star 
center and turns into zero on the core. The boundary 
condition on the star surface, r = R, is dictated by 
symmetry of the general toroidal solution of the vector 
Laplace equation which then is tested to reproduce the 
moment of inertia of a rigidly rotating solid star [l5l [27| . 
The above boundary conditions lead to the coupled alge- 
braic equations 



AtRi 



B,R- 



= 0, A e R e + BtR- 1 - 1 = R (41) 



whose solutions are 

At =Afi, Be = -Me R 2 r l+ \ Me 



R 



e+2 



R. 



2e+ 



r(42) 



In spherical polar coordinates, the nodeless toroidal field 
has only one non-zero azimuthal component 



a<t> 



0, a g = 0, 
Mr* 



(43) 



(i-C 2 ) 1/2 



dPe{Q 

dc • 



The snapshot of material node-free displacements in the 
crust undergoing torsional oscillations against immobile 
core of paramagnetic neutron star under consideration 
is pictured in Fig. 2 for quadrupole, 1 = 2, and octupole 
I = 3, overtones of this axial mode. The adopted first 



boundary condition (|37| implying that all stresses (elas- 
tic, magnetic and viscous) vanish on the core-crust inter- 
face suggests that quake-induce perturbation sets in the 
node-free torsional motions only a finite-depth crustal 
region, whereas central undisturbed region of the star 
remains at rest. In the next section the result of an- 
alytic computations are presented in the form showing 
that spectral formulas for toroidal modes entrapped in 
the crust are reduced to the above presented ones, equa- 
tions (j2~4"|) and (|2"B1 . for the global oscillations (in the en- 
tire volume of the star) when core radius tends to zero; 
this fact is regarded as a test justifying mathematical 
correctness of presented computations. 



III. SPECTRAL FORMULAE FOR THE 
FREQUENCY AND LIFETIME 

The computation of integrals defining mass parame- 
ter Ai, parameter of vibrational rigidity K and viscous 
friction T> [which has been presented in some details 
elsewhere [lj, [27j are quite lengthy but straightforward 
and, therefore, are not presented here. The mass param- 
eter can be conveniently represented in the form 
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M=4ir pR 5 

m{£) = (1 - 
x [1 - (2£-\ 
(21 +1) 2 
+ 21-1 " 

\=*±=l-h 
R 

AR = R — R r , 



£{e + i) 



3)A 2 ^ +1 



f 1)(2£- 

l\-2 



3) 



m(£), 



(44) 
(45) 



2£+3 



2^ + 3 



21 
h = 



1 

AR 

IT 7 

< A < 1. 



^2(2£+l) 



(46) 
(47) 



The A-terms in the above and foregoing equations emerge 
as result of integration along the radial coordinate from 
radius of the core-crust interface r — R c to the star ra- 
dius, r — R. The integral coefficient of viscous friction is 
given by 



V = iirr/R 
d{t) 



3^ 



1) 



(1 



21- 

A 2f+ 



d{£), 



(1 + 2) 

(e-iy 



A 2W 



(48) 
(49) 



For the lifetime we obtain 

2t„ m(£) 



-(ote) = 



R 2 



(2i + 3)(e-l) d(£) 



= ^-(50) 
P 

In Fig. 3, the fractional lifetime is plotted as a function 
of multipole degree i with indicated values of fractional 
depths h = AR/R. It shows that the higher £ the shorter 
lifetime. It is easy to see that in the limit, A = (R c /R) — > 
0, we regain the spectral formula for lifetime of global 
torsional nodeless vibrations of solid star(7l. f]~3lj 
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Alfven toroidal mode 



FIG. 3. (Color online) The fractional lifetime of torsion node- 
less oscillations of the neutron star crust damped by force of 
viscous shear stresses as a function of multipole degree I com- 
puted at indicated values of the fractional depth h of periph- 
eral seismogenic layer. 



in which by f v is understood, in this latter case, the av- 
erage kinematic viscosity of the star matter as a whole; 
the extensive discussion of this transport coefficient can 
be found in[37j. For the node-free torsional oscillations of 
solid star, the last equation has one and the same physical 
significance as the well-known Lamb formula does for the 
time of viscous damping of spheroidal node-free vibra- 
tions which in the context of neutron star pulsations has 
been extensively discussed in[38| . Regarding the problem 
under consideration we cannot see, however, how the ob- 
tained formulae can be applied to observational data on 
QPOs in SGRs. Nonetheless, their practical usefulness is 
that they can be utilized in the study of a more wide class 
of solid celestial objects such as Earth-like planets (l8l.lT9j 
and white dwarf stars [39]. 

From above it is clear that the integral coefficient of 
elastic rigidity /C e of torsional vibrations has analytic 
form similar to that for coefficient of viscous friction T>, 
namely 



p 2£ + 1 

fc e W = (l-A 2£+1 )- 1 



1 - 



(* + 2) 
{I -I) 



21+1 



(52) 
(53) 



The frequency as a function of multipole degree £ of node- 
free elastic vibrations in question v e {£) (measured in Hz 
and related to angular frequency as uj e {£) = 2iri> e (£) = 
JC e /M) is given by 



vi{t) = vi[{2i+m-i)] 



k e {£) 



(54) 



uj e = 2-KV e = — , ct = \ [^-, A = 1 — h, h = . (55) 

R V p' ' R v ' 

It is easy to see that in the limit A — > 0, we regain spectral 
formula for the frequency of global torsional oscillations, 
having the form of equation (|54p with k e (£) — m{£) = 
1. Understandably that in this latter case all material 
characteristics belong to the star as a whole. 
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FIG. 4. (Color online) Fractional frequency and period 
of nodeless torsional magneto-solid-mechanical oscillations, 
toroidal Alfven mode - oa\, entrapped in the neutron star 
crust as functions of multipole degree I computed at indi- 
cated values of the fractional depth h of peripheral seismo- 
genic layer. The value h = 1 corresponds to global torsional 
oscillations excited in the entire volume of the star. Here 
va = uja/^tt, where uia = va/R with va ~ B/^/Anp being 
the velocity of Alfven wave in crustal matter of density p and 
Pa = 2ty/loa- 



The magneto-mechanical stiffness of Alfven vibrations 
)C m can conveniently be written as 



K m = B 2 R 3 



£(£ 2 -!){£+ I) 
{2£+l)(2£- 1) 

a) = (l-A 2W )-Ml 



k m (£) 



{P - l)(2^ + 3) 



(56) 
x (57) 



1 - -£(£ + 2)(2£- l)X 2e+1 



This leads to the following two-parametric spectral for- 
mula 



vim = < 



(£ 2 1) 



21 + 3 
2£- 1 



k m 00 

m{£) ' 



UJ A = 2-kv a = — , 
ri 



D 



V A = 



(58) 
(59) 



In the forward asteroseismic analysis of QPO data rely- 
ing on this latter spectral formula, the Alfven frequency, 
va and the fractional depth of seismogenic zone, h, are 
regarded as free parameters which are adjusted so as 
to reproduce general trends in the observed QPOs fre- 
quencies. In Fig. 4, the fractional frequencies and periods 
of this toroidal Alfven mode as functions of multipole 
degree £ are plotted with indicated values of fractional 
depth of the seismogenic layer h. Remarkably, the lowest 
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overtone of global oscillations is of quadrupole degree, 
1 = 2, whereas for vibrations locked in the crust, the 
lowest overtone is of dipole degree, £ = 1, as is clearly 
seen in Fig. 5. This suggests that dipole vibration can be 
thought of as Goldstone's soft mode whose most conspic- 
uous property is that the mode disappears (the frequency 
tends to zero) when key parameter regulating the depth 

0"F spismncrpnir 7nnp \ — ^ 
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FIG. 5. (Color online) Fractional frequency of nodeless tor- 
sional Alfven oscillations of indicated overtones £ as a function 
of the fractional depth h of peripheral seismogenic layer. The 
vanishing of dipole overtone in the limit of h — ¥ 1, the case 
when entire mass of neutron star sets in torsional oscillations, 
suggests that dipole vibration possesses property typical to 
Goldstone's soft modes . 



IV. APPLICATION TO QPOS IN THE 
OUTBURST X-RAY FLUX FROM SGR 1806-20 
AND SGR 1900+14 

The basic physics underlying current understanding of 
interconnection between quasi-periodic oscillations of de- 
tected electromagnetic flux and vibrations of neutron star 
has been recognized long ago[23l l2oT |. Owing to the ef- 
fect of strong flow-field coupling, which is central to the 
propagation of Alfven waves, the quake induced pertur- 
bation excites coupled vibrations of perfectly conduct- 
ing solid-state plasma of the crust (as well as gaseous 
plasma of magnetar corona expelled from the surface by 
outburst) and frozen-in lines of magnetic field. Outside 
the star the vibrations of magnetic field lines are cou- 
pled with oscillations of gas-dust plasma expelled from 
the star surface by quake. And it is these fluctuations of 
outer lines of magnetic field, operating like transmitters 
of beams of charged particles producing coherent (curva- 
ture and/or synchrotron) high-energy radiation, are de- 
tected as QPOs of light curves of the SGRs giant flares. 

In applying the obtained spectral formulae to the fre- 
quencies of detected QPOs we examine two scenarios, 
namely, when quake-induced torsional vibrations are re- 
stored by joint action of Lorentz magnetic and Hooke's 
elastic forces and when oscillations are of pure Alfven's 



nature, that is, produced by torsional seismic vibra- 
tions of crust against core under the action of solely one 
Lorentz force of magnetic field stresses. 



A. Crust vibrations driven by combined action of 
Lorentz magnetic and Hooke's elastic forces 

In this case, the asteroseismic analysis of detected 
QPOs rests on the three-parametric spectral formula 

v 2 (£) [u A ,v e ,h] = v 2 m {£) \v A ,h] + v 2 e (£)[v e ,h]. (60) 

The suggested theoretical £-pole specification of the de- 
tected frequencies is presented in Fig. 6 for SGR 1900+14 
and in Fig. 7 and Fig. 8, exhibiting remarkable correla- 
tion between depth of seismogenic zone and fundamen- 
tal frequencies of magnetic and elastic oscillations - the 
larger AR, the higher basic frequencies of Alfvenic va 
and elastic v e vibrations. It is seen from computations 
for SGR 1806-20, that reasonable fit of data can be at- 
tained with h = 0.2 (for the star model with radius 20 
km, AR=2 km) and with h = 0.4 (AR = 5 km). 
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FIG. 6. (Color online) Theoretical fit of the QPOs frequency 
in the X-ray flux from SGR 1900+14 on the basis of three- 
parametric theoretical spectrum of frequency of torsional seis- 
mic vibrations in the crustal region of indicated fractional 
depth. 

It is worth emphasizing that at above values of h, 
the obtained here tree-parametric spectral formula much 
better match the data as compared to that for global, 
in the entire volume, vibrations studied in[16|. On this 
ground we conclude, if the detected QPOs are produced 
by seismic vibrations of peripheral region of the star un- 
der coherent action of Lorentz and Hooke's forces, then 
the depth of seismogenic layer AR should be quite large, 
somewhere in the range 0.2i? < AR < OAR. 
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FIG. 7. (Color online) The same as Fig.6, but for SGRs 1806- 
20 with h=0.2. 
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FIG. 8. (Color online) The same as Fig.6, but for SGRs 1806- 
20 with h=0.4. 



B. Lorentz-force-driven vibrations of crustal 
solid-state plasma 

It seems appropriate to note, that pure Alfven oscilla- 
tions of crustal electron- nuclear solid-state plasma about 
axis of magnetic field frozen in the immobile core have 
been studied some time ago [11] in the context of search- 
ing for fingerprints of post-glitch vibrational behavior of 
radio pulsars. In the problem under consideration, one 
can use one and the same spectral formula for the £-pole 
specification of detected QPOs which in above notations 



is written as 

v 2 (£) = vi(£)[v A ,h]. (61) 
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FIG. 9. (Color online) Theoretical description (lines) of de- 
tected QPO frequencies (symbols) in the X-ray flux during 
the flare of SGRs 1806-20 and SGR 1900+14 as overtones of 
pure Alfven torsional nodeless oscillations of crustal magneto- 
active plasma under the action of solely Lorentz restoring 
force. 

The results presented in Fig. 9 and Fig. 10 show that at 
indicated input parameters, i.e., the Alfven frequency va 
and the fractional depth of seismogenic layer h = AR/R, 
the model too adequately reproduces general trends in 
the data with fairly reasonable ^-pole specification of 
overtones pointed out by integer numbers along x-axis. 
It is seen that the low-frequency QPOs in data for SGR 
1806-20, are interpreted as dipole and quadrupole over- 
tones: ^(o»i) = 18 and ^(oa|) = 26 Hz. And the high- 
frequency kilohertz vibrations with 627 Hz and 1870 Hz 
are unambiguously specified as high-multipole overtones: 
j/(oa^ =42 ) = 627 Hz and ^(oa 1 ^^) = 1870 Hz. However, 
in this latter scenario of Lorentz-force-dominated vibra- 
tions the best fit of data is attained at fairly large value of 
fractional depth, h — 0.6, which is much larger than the 
expected depth of the crust. In our opinion, this result 
may be regarded as indication to that the detected QPOs 
are formed by coherent vibrations of crustal solid-state 
plasma and plasma of magnetar corona. 

V. CONCLUDING REMARKS 

Ever since identification of pulsars with rapidly rotat- 
ing neutron stars it has been argued [13, EJ that two 
key properties of these compact objects - (i) the degen- 
eracy of neutron (non-conducting) Fermi-matter whose 
pressure opposes the pressure of self-gravity, and (ii) a 
highly stable to decay super-strong magnetic fields - can 
be reconciled, if poorly conducting neutron-dominated 
stellar matter, constituting the neutron star cores, has 



10 



5000 




/ 



FIG. 10. (Color online) Same as Fig.9 but for SGRs 1806-20. 



been brought to gravitational equilibrium in the per- 
manently magnetized state. The most plausible is the 
state of Pauli's paramagnetic saturation with spin mag- 
netic moments of neutrons polarized along the axis of fos- 
sil field inherited from massive progenitor and amplified 
in magnetic-flux-conserving core- collapse supernova^. 
This idea is central to the considered two-component, 
core-crust, model of paramagnetic neutron star whose 
less dense and highly conducting, metal-like, material 
of the crust is considered as a solid-state, electron- 
nuclear, plasma pervaded by frozen in the core magnetic 
field. This difference between electrodynamic properties 
of core and crust matter (permanently magnetized non- 
conducting core and perfectly conducting non- magnetic 
crust) suggests that magnetic cohesion between massive 
core (permanent magnet) and crust (metal-like material) 
should plays central role in seismic activity of the star. 
Working from such an understanding, we have computed 
frequency spectra of node-free torsional oscillations of 
crust against immobile core under the action Lorentz 
and Hookc restoring forces and damped by Newtonian 
viscous force. As a trial function of oscillating mate- 



rial displacements we have used the node-free toroidal 
field computed from vector Laplace equation. The ob- 
tained spectral formulas are of some interest in their own 
right because they can be applied to more wide class of 
celestial objects. In this work we applied the obtained 
analytic frequency spectra to magnetars, highly magne- 
tized quaking neutron stars whose bursting seismic activ- 
ity is commonly associated with release of magnetic field 
stresses. Focus was laid on forward asteroseismic anal- 
ysis of fast X-ray flux oscillations during the giant flare 
of SGR 1900+14 and SGR 1806-20 and, thus, assuming 
that these oscillations are produced by torsion vibrations 
of crustal solid-state plasma about axis of dipole mag- 
netic field frozen in the immobile permanently magne- 
tized core. In so doing we have investigated two cases 
of post-quake vibrational relaxation of the star, depend- 
ing on restoring forces. In first case, the analysis of data 
has been based on assumption that detected QPOs owe 
their existence to node-free torsional vibrations of crust 
against core restored by joint action of Lorentz magnetic 
and Hooke's elastic forces. And we found that obtained 
three-parametric spectral formula provides much better 
fit of data than two-parametric frequency spectrum of 
global vibrations (Bastrukov et al 2009a). The consid- 
ered second scenario presumes that vibrations are domi- 
nated by solely Lorentz restoring force of magnetic field 
stresses. We found that obtained two-parametric fre- 
quency spectrum can too be fairly reasonably reconciled 
with detected QPOs frequencies. All the above lead us 
to conclude that Lorentz restoring force of magnetic field 
stresses plays decisive part in quake-induced torsional vi- 
brations of crustal solid-state plasma of magnetars. 

The authors are grateful to Dima Podgainy (JINR, 
Dubna) for helpful assistance and referee and reviewer 
for suggestions clarifying the subject and result of pre- 
sented investigation. 
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